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ABSTRACT 

Transition disc systems are young stars that appear to be on the verge of dispersing 
their protoplanetary discs. We explore the nature of these systems by comparing the 
stellar accretion rates M„ and disc masses of transition discs and normal T Tauri 
stars in Taurus and Ophiuchus. After controlling for the known dependencies of M, 
and Md on age, M* on stellar mass, and Md on the presence of stellar or sub-stellar 
companions, we find that the normal T Tauri stars show a trend of M* increasing 
with Md. The transition discs tend to have higher average disc masses than normal T 
Tauri stars as well as lower accretion rates than normal T Tauri stars of the same disc 
mass. These results are most consistent with the interpretation that the transition 
discs have formed objects massive enough to alter the accretion flow, i.e., single or 
multiple giant planets. Several Ophiuchus T Tauri stars that are not known transition 
disc systems also have very low accretion rates for their disc masses. We speculate on 
the possible nature of these sources. 

Key words: planets and satellites: formation - protoplanetary discs - stars: forma¬ 
tion - 


1 INTRODUCTION 

The population of young stellar objects known as transition 
discs have unusual spectral energy distributions (SEDs) that 
suggest that some portion of the inner disc has become op¬ 
tically thin (Strom et al. 1989; Skrutskie et al. 1990). The 
radial distribution of the dust disc has been inferred from de¬ 
tailed modeling of the SEDs, which have significant excesses 
beyond 10pm, like those of normal T Tauri stars, but weak 
excesses at mid-infared wavelengths (e.g., D’Alessio et al. 
2005; Calvet et al. 2005). The optically thin region inferred 
from the modeling—a central “hole” in the dust distribu¬ 
tion or an annular “gap” within a radius ry of 15 AU to 
> 50 AU depending on the source—has been confirmed by 
submillimeter interferometric imaging of bright, nearby sys¬ 
tems (Andrews et al. 2011). These properties suggest that 
transition discs are in the process of dispersing their discs 
(see Espaillat et al. 2014 for a review). 

Several scenarios have been proposed to explain the 
properties of transition discs. An early suggestion was that 
grain growth and planetesimal formation, the first steps to¬ 
ward planet formation via core accretion, might reduce the 
optical depth of the inner disc by reducing the dust cross- 
sectional area (Strom et al. 1989; see also Dullemond & Do- 
minik 2005; Tanaka et al. 2005). The formation of giant 
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planets is also expected to carve a gap or inner cavity in the 
disc, rendering some or all of the inner disc optically thin 
(Skrutskie et al. 1990; Marsh & Mahoney 1992, Lubow et al. 
1999; Calvet et al. 2002; Rice et al. 2003; Quillen et al. 2004, 
Calvet et al. 2005, D’Alessio et al. 2005; Rice et al. 2006). 
Alternatively, a photoevaporative wind, launched by the ir¬ 
radiation of the disc surface by stellar UV and X-rays, may 
also deplete the inner disc, by limiting the ability of material 
accreting from the outer disc to resupply the material in the 
inner disc (Clarke et al. 2001; Alexander et al. 2006; Alexan¬ 
der & Armitage 2007; Owen et al. 2010, 2011; Gorti et al. 
2009). These interpretations are not mutually exclusive. For 
example, photoevaporation may disperse a disc more rapidly 
once a gap forms as a result of giant planet formation (e.g., 
Rosotti et al. 2013). Transition discs may also be produced 
through multiple pathways, with different processes poten¬ 
tially dominating on different timescales (e.g., Williams & 
Cieza 2011). 

Because multiple processes can produce transition disc 
SEDs, other information (beyond SEDs) has been used in 
efforts to distinguish among the potential explanations. De¬ 
tailed studies of individual transition discs have reported 
evidence for an orbiting companion (e.g., Kraus & Ireland 
2012; Huelamo et al. 2011) or its circumplanetary disc (Brit¬ 
tain et al. 2013, 2014). Other studies have probed the spatial 
and velocity distribution of disc gas to look for evidence of 
photoevaporative outflows (e.g., Pascucci & Sterzik 2009; 


Sacco et al. 2012) or the radial truncation of the gaseous 
disc (e.g., Najita et al. 2008; Dutrey et al. 2008). The latter 
is not expected in the planetesimal scenario. 

An alternative way to distingnish among the possi¬ 
ble scenarios is to compare the demographic properties of 
populations of transition discs with those of normal (non¬ 
transition) T Tauri stars. Whereas recent studies have em¬ 
ployed stellar accretion rates alone or in combination with 
properties such as X-ray luminosity and stellar mass for this 
purpose (Kim et al. 2009, 2013; Espaillat et al. 2012; Fang 
et al. 2013; Manara et al. 2014), the combination of stellar 
accretion rates and disc masses has also been suggested as a 
potentially powerful way to distinguish among the scenarios 
described above (Najita et al. 2007; Alexander & Armitage 
2007; see also Cieza et al. 2008, 2010, 2012; Mendigutfa et 
al. 2012). To summarize the basic idea: 

(a) In the grain and planetesimal growth scenario, the dust 
component of the disc within is altered but not the gas; 
the stellar accretion rate is therefore expected to be similar 
to that of normal T Tauri stars. Because planet formation is 
a common outcome of disc evolution (e.g., Fressin et al. 2013; 
Dong & Zhu 2013; Petigura et al. 2013), and planetesimal 
formation is the first step toward planet formation via core 
accretion, planetesimals are expected to form in discs of all 
masses. 

(b) In contrast, when a forming Jovian-mass planet clears 
a gap in the disc, most of the accretion flow from the outer 
disc accretes onto the planet (in streams), and the stellar 
accretion rate is reduced to ~ 0.1 of its original T Tauri 
value (Lubow et al. 1999; but see also Zhu et al. 2011). 
Because a Jovian-mass planet represents a significant mass 
reservoir relative to the median T Tauri disc mass (Andrews 
et al. 2013), planets massive enough to open a gap (e.g., > 
IMj) are expected to form in discs with higher than average 
masses. 

(c) The accreting planet could eventually reach a mass high 
enough (-^ 5Mj) to suppress flow past the planet. When 
the remaining inner disc material accretes on to the star, 
an inner cavity in the gas and dust distribution of the disc 
would remain, with no further stellar accretion (e.g., Lubow 
et al. 1999; Lubow & D’Angelo 2006). Forming such high 
mass planets is favored in higher mass discs. 

(d) If a photoevaporative wind can remove material from the 
disc faster than it can be replenished by accretion, the inner 
disc is starved of gas and dust, and stellar accretion drops 
to a low or negligible rate. This situation is more readily 
achieved in low mass discs (all other factors being equal; 
e.g., stellar irradiation fields) with initially modest accretion 
rates (e.g., Clarke et al. 2001; Alexander & Armitage 2007; 
Owen et al. 2010). 

A previous study of stellar accretion rates and disc 
masses found evidence for a modest trend between stellar 
accretion rate and disc mass among single, normal T Tauri 
stars in Taurus (Najita et al. 2007). In addition, the Taurus 
transition discs were found to have disc masses 4 times 
larger than the normal T Tauri stars and stellar accretion 
rates ^ 10 times lower than those of normal T Tauri stars 
of the same disc mass. These properties are most consistent 
with the interpretation that the Taurus transition discs are 
forming Jovian mass planets. 

Here we extend the Taurus study by examining the stel¬ 
lar accretion rates M* and disc masses Md of transition discs 


and T Tauri stars in Ophiuchus in combination with the ear¬ 
lier Taurus results. Ophiuchus is one of the few star forming 
regions in which disc masses are available for a significant 
number of young stellar objects. Our study differs in several 
ways from the approach taken by Najita et al. (2007). The 
transition disc sample studied here is restricted to sources 
with evidence for large inner holes, i.e., sources with either 
a large MIR dip in the SED or a resolved cavity in submil¬ 
limeter imaging of the source. Because of the large range 
of stellar spectral types in the Ophiuchus sample, we also 
correct for the dependence of stellar accretion rate on stel¬ 
lar mass. Our data sets are described in §2 and our analysis 
and results in §3. In §4, we discuss the location of transition 
discs in the M«-Md plane and compare our results to those 
of previous studies. Our conclusions are summarized in §5. 


2 DISK MASSES AND STELLAR ACCRETION 
RATES 


For both the Taurus (Tau) and Ophiuchus (Oph) samples, 
we use stellar accretion rates and disc masses from the liter¬ 
ature, with a few emendations as described below. Because 
we intend to compare the results for Oph with those previ¬ 
ously obtained for Tau, the disc masses and stellar accretion 
rates for the two samples are treated as similarly as possible. 

Millimeter continuum fluxes have been reported previ¬ 
ously for many Oph sources. To incorporate these measure¬ 
ments in a homogeneous way in deriving disc masses, we 
follow the approach of Andrews et al. (2013) in their anal¬ 
ysis of Tau disc masses. For each disc, we compiled and fit, 
using a power law model, all the available continuum mea¬ 
surements in the literature at wavelengths from LOOqm to 
3 mm (Andre & Montmerle 1994; Dent et al. 1998; Nuern- 
berger et al. 1998; Motte et al. 1998; Stanke et al. 2006; 
Andrews & Williams 2007a, b; Patience et al. 2008; An¬ 
drews et al. 2009, 2010; Ricci et al. 2010). These model fits 
and their associated uncertainties were extrapolated to es¬ 
timate a flux density at a fiducial wavelength of 1.33 mm. If 
there was only 1 datapoint available in this range, we ap¬ 
plied a mean power-law scaling to estimate the 1.33 mm flux 
density, where oc 

We then calculated disc masses assuming the standard 
optically thin, isothermal approach, i.e.. 


Afd.thin 


KuB,y{Td) ’ 


( 1 ) 


for a fixed distance d (125 pc for Oph, 140 pc for Tau) 
and opacity (Ki. 33 mm(dust) = 2.3cm^/g; dust-to-gas ratio 
= 0.01). In these calculations, the average dust temperature 
for each source was computed from a simple estimate of the 
incident irradiation field, such that Td = 25K(L«/L 0 )^/^. 
This approximation assumes that the dust in discs is located 
at a common effective orbital distance, which is a reasonable 
assumption when spatially resolved information is unavail¬ 
able (see Andrews et al. 2013 for details). 

For Tau, we adopted the disc masses tabulated by An¬ 
drews et al. (2013) using the above method. For Oph, we 
determined Ma.thin using stellar luminosities from McClure 
et al. (2010). For both the Tau and Oph samples, we also 
investigated the impact of adopting Md values estimated 
from fitting the SED and resolved millimeter imaging data 



Afd,flt- These are available for the brighter Tau and Oph 
sources (Andrews & Williams 2007b; Andrews et al. 2009, 
2010 , 2011 ). 

The primary source of Oph stellar accretion rates is the 
large sample reported by Natta et al. (2006). Because of 
the high extinction toward pOph, Natta et al. determined 
accretion rates from the IR hydrogen lines (Pa/3 or Bry) 
using the calibration of Muzerolle et al. (1998) to convert 
the IR line luminosities to accretion luminosities. For con¬ 
sistency with the adopted Oph disc masses, which involve 
a correction for L*, we rederived the stellar accretion rate 
for each Oph source from the Pa /3 equivalent width and J- 
band magnitude reported by Natta et al. (2006) using the 
extinction and L* from McClure et al. (2010). The stellar 
masses used in the accretion rate determination were derived 
from the Siess et al. (2000) stellar evolutionary models us¬ 
ing the approach and assumptions adopted by Andrews et 
al. (2013). 

A few additional Oph accretion rates were adopted from 
Espaillat et al. (2010), Eisner et al. (2005), and Valenti et 
al. (1993). The accretion rate determined by Espaillat et 
al. (2010) for DoAr44 (also known as ROXs-44) was de¬ 
rived from the 17-band excess reported by Bouvier & Ap- 
penzeller (1992) following the Gullbring et al. (1998) method 
(Andrews et al. 2011; C. Espaillat 2013, personal commu¬ 
nication). Eisner et al. (2005) report an accretion rate for 
WaOph6 (V2508 Oph) that is based on a 17-band excess de¬ 
rived from UBVRI photometry and high resolution optical 
spectroscopy following the method of Gullbring et al. (1998). 
Because the Muzerolle et al. (1998) relations used by Natta 
et al. (2006) were originally calibrated with accretion rates 
derived via the Gullbring method, the set of M* values used 
here should be on a consistent footing. The properties of the 
Oph sample are shown in Table 1. 

For the Tau sample, we adopted stellar accretion rates 
from Najita et al. (2007) and also included a few additional 
stellar accretion rates from Valenti et al. (1993; GI Tau, 
HN Tau) and Hartigan et al. (1995; HK Tau). As in Najita 
et al. (2007), the values from Hartigan et al. (1995) were 
scaled downward by a factor of 0.17 for consistency with 
the stellar accretion rates of Hartmann et al. (1998)Q In 
order to investigate the impact of the known trend between 
M* and M*, we also adopted the stellar masses reported by 
Andrews et al. (2013) using the Siess et al. (2000) models. 
The properties of the single star Tau sample are shown in 
Table 2. 

The Oph sample includes both sources in the ~ 2 Myr 
old surface population and the surrounding region (Wilking 
et al. 2005) as well as the younger embedded sources in the 
Oph core (< 1 Myr; Luhman & Rieke 1999). In comparison, 
Tau has a median age of ~ 2 Myr with a distribution between 
1 Myr and a few Myr (Kenyon & Hartmann 1995; Hartmann 
2001; Luhman et al. 2003; Andrews et al. 2013). 


3 RESULTS 

Despite the large samples of Oph stellar accretion rates and 
disc mass estimates reported by Natta et al. (2006) and 

^ In Najita et al. (2007) we also scaled the stellar accretion rates 
of White & Ghez (2001) up by 2.4. 
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Figure 1. Stellar accretion rates (in MQyr“^) and disc masses 
(in Mq) for Ophiuchus sources. Arrows indicate upper limits. 
Closed (open) circles designate sources without (with) known 
companions within 1". Binaries have lower disc masses than sin¬ 
gle stars. Sources with spectral types earlier than MO (symbols 
with a red rim) have higher average stellar accretion rates than 
sources with later spectral types. 


Andrews & Williams (2007a) (~140 sources in each) the 
overlap between the two samples is limited (~ 40 sources). 
The limited overlap is partly due to different areal coverage. 
Natta et al. (2006) studied the Oph core, whereas Andrews 
& Williams (2007a) studied a larger area. In addition, both 
samples include a signihcant fraction of sources that are ei¬ 
ther Glass HI or Glass I (25%-35%). Figure 1 shows sources 
that have detections or upper limits (indicated by arrows) 
for both quantities. Sources with a known stellar companion 
within 1” (Reipurth & Zinnecker 1993; Barsony et al. 2005; 
Ratzka et al. 2005) are shown as open circles. Sources with¬ 
out such a companion are regarded as single stars (solid sym¬ 
bols). The binaries have lower average disc masses than the 
single stars, consistent with similar trends noted elsewhere 
in the literature (e.g., Jensen et al. 1994, 1996; Andrews & 
Williams 2005; Harris et al. 2012). The presence of a stel¬ 
lar companion is believed to dynamically alter the structure 
of the disc, reducing its mass. We therefore excluded the 
binaries from further consideration. 

In Figure 1, sources with later spectral types (later than 
Ml) have lower average M* than earlier types, a result that 
is not surprising given the well-known trend of increasing 
M* with M* (e.g., Muzerolle et al. 2003; Galvet et al. 2004; 
Natta et al. 2006; Fang et al. 2009; Rigliaco et al. 2011; 
Manara et al. 2012). This trend was not an issue in the 
earlier Taurus study, because the Tau sample spans a more 
limited range of spectral types (K2-M2; ~ 80% are K3-M0 
or 0.5 — 1.5Mq). Because the trend of M* with M, has 
a large scatter at 1 Myr age (~ 2 orders of magnitude 
at a given M*; e.g.. Fang et al. 2013), the trend is only 
apparent over a mass range of > 1 dex (e.g., Muzerolle et al. 
2003; Manara et al. 2012), as in our Oph sample (A1-M6; 
-O.O7-4M0). 

We therefore corrected the Oph sample for the trend 



between M* and M* in order to include its larger range of 
spectral types and maximize its sample size. The derived 
Oph accretion rates M«o were rescaled as log M*,norm = 
logM*o — 1.3log(Mt/O.7M0) using the slope from Manara 
et al. (2012) for a 1 Myr population and our estimated stel¬ 
lar masses (section 2). In exploring alternative approaches, 
we found that the effect of the rescaling on our results is 
insensitive to the details of which stellar evolutionary tracks 
are used in estimating stellar masses or whether the appro¬ 
priate relation has a steeper slope (e.g., M* oc M* as in 
Muzerolle et al. 2003). We also explored but did not adopt 
a correction for a possible trend of Ma oc M* (Andrews et 
al. 2013), because it did not affect our results. This is under¬ 
standable because Ma and M* are not strongly correlated 
over the primary mass range of our sample (logM*= -0.65 
to 0.3; Andrews et al 2013). 

Figure 2a shows the rescaled Oph accretion rates plot¬ 
ted against Ma,thin (black dots) where, for clarity, sources 
with upper limits on their disc mass are not shown. Sources 
with a transition disc SED and/or a submillimeter cavity 
(Espaillat et al. 2010, McClure et al. 2010; Andrews et al. 
2011) are indicated by blue ‘-I-’ symbols. Single T Tauri stars 
(gray dots) and transition discs (blue ‘x’ symbols) in Taurus 
are also shown (see Najita et al. 2007 and section 2), where 
we have applied the same rescaling of M* as a function of 
M* that was used for the Oph sample. The Oph transition 
discs fall between the regions occupied by the Tau transition 
discs and the normal T Tauri stars. 

Most of the sources marked as transition discs in Ta¬ 
ble 1 and Figure 2 have prominent mid-infrared (MIR) dips 
in their SEDs. Because there are only 2 such Oph objects 
with M* measurements (SR21A, DoAr44; McClure et al. 
2010), we also included Oph sources that have a submil¬ 
limeter cavity (but lack a MIR SED dip) for which M* 
is available (SR24S, WSB60; Andrews et al. 2011). These 
sources (both with and without MIR SED dips) have sub¬ 
millimeter cavity radii of 15-36 AU (Andrews et al. 2011). 
Because submillimeter imaging studies have been limited to 
the brighter sources (> 75mJy at SSOpm or Ma > 0.005Mq 
for the Oph SMA sample; Andrews et al. 2010), the sample 
of cavity sources may be artificially biased to higher disk 
masses; i.e., lower mass disks may possess cavities that have 
not yet been identified. In contrast, selecting transition discs 
by their prominent MIR SED dip should not introduce a 
similar bias. For example, a system like CoKu Tau/4, which 
has a prominent MIR SED dip, has a very low disk mass 
(- 0.001 M©; D’Alessio et al. 2005). 

The Tau transition disc sample differs from that of Na¬ 
jita et al. (2007), who defined transition discs more broadly 
based on weak or no excess shortward of lOpm and a sig¬ 
nificant excess at longer wavelengths. Because such a broad 
definition is debated, we have adopted a more conservative 
definition here by selecting only objects with prominent MIR 
dips in their SEDs (DM Tau, GM Aur, UX Tau, LkCal5; 
Calvet et al. 2005; Espaillat et al. 2010). The remaining 
sources in the Tau sample are classified as normal T Tauri 
stars. 

To explore the impact of our simple prescription for disc 
mass (eq. 1), Figure 2b shows the same information except 
that the Tau and Oph disc masses estimated from radiative 
transfer models that simultaneously fit SEDs and resolved 
millimeter-wavelength visibilities are used when avail- 
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Figure 2. Top (a): Stellar accretion rates (in MQ,yr~^) and disc 
masses (in Mq) of T Tauri stars in Ophiuchus (black dots) and 
Taurus (gray dots) compared with transition discs in Ophiuchus 
(blue ’+’) and Taurus (blue ’x’). Red circles indicate the Ophi¬ 
uchus outlier population, i.e., T Tauri stars with normal SEDs 
that have very low accretion rates for their disc masses. For clar¬ 
ity, only non-binaries and sources with detections are shown. 
Stellar accretion rates have been adjusted for the known depen¬ 
dence of M* on M*. Arrows indicate sources with upper limits 
on their stellar accretion rates. The dotted line is a linear fit to 
the normal T Tauri stars; the dashed line, located 1-cr below the 
dotted line, divides the plane into two regions (see text for de¬ 
tails). Bottom (b): Similar to the top panel except that disc mass 
estimates from fits to resolved millimeter-wavelength visibilities 
Md. fjt are used when available. In either plot, transition discs 
have lower than average accretion rates for their disc mass. 


able (Andrews & Williams 2007b; Andrews et al. 2009, 2010, 
2011). Of the 13 Oph sources with Mdfit values, roughly half 
have values within 60% of the Mayhin value. The other 
half have Md.at values ~ 5 times larger. The large difference 
for sources such as WaOph6 and DoAr25 occurs because the 
dust associated with these sources is colder and the disc is 




larger in size than one would infer from the simple optically 
thin scaling used for Md,thin- A similar result is found for 
the Tau sources. 

Note that the Tau accretion rates and disc masses 
shown in Figure 2a differ from those adopted in Najita et 
al. (2007). Here the accretion rates are scaled by M«, and 
the optically thin disc masses assume a common effective 
orbital distance rather than a common average tempera¬ 
ture for the emitting dust grains. That is, the characteristic 
(mass-weighted) temperature of the emitting grains, when 
averaged over the radial structure of the disc, depends on 
L«; in contrast, the Andrews & Williams (2005) disc masses 
used in Najita et al. (2007) assumed a characteristic grain 
temperature of 20 K for all sources (for details see section 
3.2.2 of Andrews et al. 2013). 

As shown in Figure 2, with either set of disc mass es¬ 
timates, the normal (single) T Tauri stars in Oph overlap 
the region of the diagram occupied by the normal (single) 
T Tauri stars in Tau. Among the sources with measured 
stellar accretion rates and disc masses, there is a trend of 
M» increasing with Md, albeit with significant dispersion. A 
similar result was noted previously for the Tau sample alone 
(Najita et al. 2007). A trend of M* increasing with Md is 
not surprising; in a-viscosity discs (e.g., Shakura & Sunyaev 
1973) the accretion rate is proportional to the disc column 
density (or mass). 

Time variability and measurement uncertainty in the 
stellar accretion rate likely account for some of the observed 
scatter about the mean trend. Stellar accretion rates of in¬ 
dividual objects are known to vary by a factor of ~ 2 (e.g., 
Hartigan et al. 1991). Gullbring et al. (1998) reported typ¬ 
ical measurement uncertainties of 0.5 dex in their stellar 
accretion rates; Natta et al. (2006) did not report their un¬ 
certainties. 

The assumptions made in deriving disc masses are also 
a likely source of dispersion. Disks may be partially opti¬ 
cally thick or have different temperature distributions than 
assumed in the prescription of Andrews et al. (2013). In ad¬ 
dition, Md measurements are insensitive to solids that have 
grown into large grains (> cm). Dispersion in the solid size 
distribution would mean dispersion in the extent to which 
the millimeter continuum traces the gas mass. Unrecognized 
binarity (and its impact on either or both of M* and Md) is 
yet another potential source of dispersion, although more so 
for Oph than for Tau. The binarity of the latter population 
is now fairly well known (e.g., Kraus et al. 2011). 

Several additional Oph sources have low M« for their 
Md but lack obvious MIR dips in their SED (Fig. 2, red 
open circles). DoAr 25, GY284, and WL14 have upper lim¬ 
its below logM» = —9.5, and LFAM3 has the lowest Oph 
logM* measurement in the plot. WL14 and LFAM3 have 
high estimated extinction values (A„ > 15), with their stel¬ 
lar accretion rates possibly underestimated as a result. We 
discuss this issue and the possible nature of these “outlier” 
sources in the next section. 

We also find that the known Oph transition discs oc¬ 
cupy roughly the same region of the M«-Md plane as the 
Tau transition discs. That is, the Oph transition discs have 
lower accretion rates for their disc masses than do normal T 
Tauri stars in Tau and Oph. In Figure 2a, the median log M» 
for the normal T Tauri stars (outliers excluded) and transi¬ 
tion discs is -7.6 and -8.4, respectively; the median logMd 


is -2.3 and -1.7, respectively. If the outliers were included in 
the normal T Tauri population, the median log M* would 
decrease slightly to —7.8; the median log Ma would be un¬ 
changed. 

With the addition of the Oph sources, we can describe 
more quantitatively the trend between log Md and log M* for 
the normal T Tauri population (outliers excluded) in Fig¬ 
ure 2a. The Kendall r for the individual Tau and Oph and 
the combined (Tau+Oph) normal T Tauri samples is 0.34, 
0.46, and 0.38, respectively, with corresponding 2-sided p- 
values of 0.96%, 2.8%, and 0.03%. The reduced p-value of 
the combined sample illustrates the value of the Oph sources 
in improving the statistical significance of the trend. Includ¬ 
ing the outliers in the normal T Tauri population reduces 
the Kendall r to 0.23 (p-value of 2.3%), consistent with the 
visual impression of these sources as outliers. 

We estimated the slope of the trend using Brandon 
Kelly’s linmix^err idl routine (Kelly 2007), which takes a 
Bayesian approach to linear regression given errors in both 
quantities. Assuming errors of 0.2 dex in log Md,thin and 0.3 
dex in log M*,norm, we find 

log M*,norm = -5.15(±0.75) -f 1.08(±0.3) log Md.thin (2) 

with a dispersion about this trend of (T(log M»,norm) = 0.5. 
If the Tau and Oph normal T Tauri star populations are 
examined separately, they have similar slopes (-^ 1) and dis¬ 
persion (~ 0.5 dex). The M*,norm normalization constant 
for the Oph normal T Tauri population is larger by 0.75 dex 
than that of the Tau population, consistent with the visual 
impression, although the uncertainty in the values of the 
constants is large. 

The small sample size of the transition discs makes it 
difficult to compare in detail the distributions of the normal 
T Tauri stars and transition disc populations in the log Ma- 
log M« plane. For example, if we were to ht a linear trend to 
the transition disc sample for comparison with the fit for the 
normal T Tauri stars (eq. 2), the uncertainty on the slope 
and constant would be too large for a meaningful compar¬ 
ison. However, we can explore whether the distributions of 
the two populations in the M«-Md plane are grossly similar 
by dividing the plane into two regions and seeing whether 
there is a statistical difference in the distribution of the two 
populations across the two regions. Using the distribution of 
normal T Tauri stars as a guide, since it is the better char¬ 
acterized population, we divide the plane at a line located 
1 -cr below the mean log Md,thin-log M», norm relation for the 
Tau+Oph normal T Tauri stars (eq. 2), i.e., at 

log M*,norm = -5.65 + 1.08 log Md.thin (3) 

(Fig. 2a, dashed line). As expected, many more of the normal 
T Tauri stars in the combined sample lie above the line (32 
sources) than below the line (12 sources). In contrast, there 
are no transition discs above the line and 8 below. For the 
purpose of apportioning the sample, upper limits in logM* 
are treated as detections, although the results are insensitive 
to how they are treated. 

We can use a contingency table analysis (e.g., Bussmann 
et al. 2011) to see whether these values are consistent with 
the null hypothesis that the populations are similarly dis¬ 
tributed above and below the line. Fisher’s exact test, which 
is a useful approach when the sample sizes are small, gives a 
probability of 0 that the T Tauri stars and transition discs 



are similarly distributed. The result is insensitive to where 
the boundary (eq. 3) is drawn. If the boundary were a factor 
of 2, 2.5, or 3 times higher in logM*, the probability would 
still be only 0.002, 0.006, or 0.014 that the T Tauri stars and 
transition discs are similarly distributed. The boundary in 
the last case is coincident with the locus of the mean trend 
(Fig. 2a, dotted line). If the boundary were a factor of 2 
or 3 lower in logM*, the probability remains low at 0.001 
or 0.002 respectively. Including the outliers in the normal T 
Tauri population raises the above probabilities slightly (by 
a factor of ~ 2) but does not change our finding that the 
transition discs and normal T Tauri stars are not similarly 
distributed in the M*-Md plane. As a contrasting example, 
we can also compare the distributions of the Tan and Oph 
normal T Tauri stars across the boundary (eq. 3). The prob¬ 
ability that the two populations are similarly distributed is 
21 %. 


4 DISCUSSION 

As described in the previous section, when combined with 
the Tan sample, the Oph sample enhances the statistical 
signficance of the previously reported trend of M* vs. Md 
among single, normal T Tauri stars. The Oph transition 
discs, like the Tan transition discs, are found to favor the 
lower right region of the M«-Md plane, i.e., transition discs 
tend to have lower M* for a given Md and a higher Md (for 
a given M«) than normal T Tauri stars. 

Planetesimal formation does not provide a ready expla¬ 
nation for this result. Naively, we would expect that systems 
forming planetesimals would have accretion rates similar to 
those of normal T Tauri stars, because planetesimal forma¬ 
tion alters the size distribution of the disc solids, but it is 
not known to alter the properties of the disc gas and stel¬ 
lar accretion rate. While this result was obtained previously 
for the Tan transition disc sample, it also appears to ap¬ 
ply at the younger average age of Oph, when planetesimal 
formation, as the first step in planet formation in the core 
accretion scenario, would be even more likely to dominate 
over the other proposed explanations for transition discs (gi¬ 
ant planet formation, photoevaporation). The lack of transi¬ 
tion discs meeting our expectations for where planetesimal- 
forming systems would lie in the M*-Md plane, for both the 
Tau and Oph samples, argues against planetesimal forma¬ 
tion as a significant pathway to form transition discs at the 
ages of nearby star forming regions (1 — 3Myr). 

This interpretation is consistent with the results of re¬ 
cent theoretical studies that find it difficult to replicate the 
properties of transition discs through grain evolution alone. 
Early models of grain evolution that considered sticking 
without fragmentation were able to produce MIR SED dips 
like those of transition discs (Dullemond & Dominik 2005; 
Tanaka et al. 2005). However, recent, more realistic grain 
evolution models that include coagulation, fragmentation, 
radial drift, gas drag, and turbulent mixing (Birnstiel et al. 
2012) find that while discs can develop the MIR dips that 
characterize transition disc SEDs, enough non-fragmenting 
centimeter-sized particles remain behind that these systems 
would not show the submillimeter cavities that also char¬ 
acterize transition discs. A similar difficulty is found even 
under the simplified conditions considered in earlier models 


(perfect sticking in the absence of radial draft; e.g., Dulle¬ 
mond & Dominik 2005). 

The measurable M» and higher average Md of transi¬ 
tion discs compared to normal T Tauri stars also suggests 
that photoevaporation is not the dominant pathway to a 
transition disc at Taurus-Ophiuchus age. Photoevaporation 
is more likely to create an inner cavity in the disc in systems 
with low M« and low Md, if all other factors are equal (e.g., 
Alexander & Armitage 2007). 

The low M« of transition discs for their disc masses rel¬ 
ative to normal T Tauri stars is most consistent with the 
formation of objects massive enough to alter the disc ac¬ 
cretion flow, i.e., giant planet formation. Interestingly, this 
is the case for both systems with transition disc SEDs and 
those with submillimeter cavities. Sources such as DoAr44, a 
transition disc with a prominent MIR SED dip, are hypoth¬ 
esized to have formed giant planetary companion(s) (e.g., 
Espaillat et al. 2010). 

Why do we find low accretion rates for transition discs 
both with and without MIR dips in their SEDs? Perhaps 
the conditions needed to reduce the stellar accretion rate 
are less stringent than those needed to render the inner disc 
optically thin in the MIR. Removing small dust grains from 
the inner disc enough to render it optically thin while also 
allowing continued accretion onto the star is a recognized 
challenge (Rice et al. 2006; Zhu et al. 2012). Whether a sys¬ 
tem with a submillimeter cavity also develops an optically 
thin inner disc in the MIR (i.e., a transition disc SED) may 
therefore require the formation of higher mass companions 
and/or a larger number of companions (Zhu et al. 2011; 
Dodson-Robinson & Salyk 2011). For example, a more mas¬ 
sive companion or more of them may induce a larger pres¬ 
sure bump at the inner edge of the outer disc, enhancing the 
efficiency of dust filtering (Rice et al. 2006). Radiation pres¬ 
sure from a high mass accreting planet can also play a role 
in creating a transition disc SED (Owen 2014). In compar¬ 
ison, the formation of a single, low mass giant planet may 
be sufficient to create a submillimeter cavity and to impact 
the disc accretion rate (e.g., Lubow et al. 1999). 

In §3, we noted a few outliers among the Oph sources 
in the M»-Md plot, i.e., sources with normal T Tauri SEDs 
but low accretion rates for their disc mass. For some sources, 
their high Ay may have led to an underestimate of their ac¬ 
cretion rates (WL14, LFAM3). In sources with very high 
reddening, such as highly-inclined discs or protostars, the 
scattered light component can dominate the observed light 
into the near-infrared. Because the colors then appear bluer 
than would be produced by the extinction along the line of 
sight to the central source, the typical methods of determin¬ 
ing the extinction using observed near-infrared colors and 
assuming an intrinsic color from a known spectral type or 
the classical T Tauri star locus will underestimate the ex¬ 
tinction. Subsequent dereddening of the observed line flux 
will lead to an underestimate of the true line luminosity, 
and by extension, M«. Correcting for this effect for a given 
source depends on the scattering fraction, which is difficult 
to determine. 

Another possibility is that some or all of the low ac¬ 
cretion rate sources may be unrecognized transition discs in 
which giant planet formation either fij has occurred at disc 
radii too small to probe with SEDs or submillimeter imaging 
or fuj is insufficiently advanced (e.g., the planet is too low 



in mass or too few in number) to render portions of the disc 
optically thin. This scenario may be particularly relevant to 
the two outliers that have only modest extinction (DoAr25 
at Av = 3.4 and GY284 at Ay = 6). 

Andrews et al. (2008) previously noted that DoAr25 
has a low stellar accretion rate for its disc mass. Al¬ 
though DoAr25 has one of the highest measured disc masses 
(O.I4M0), stellar accretion estimates range from < 2 x 
lO~^°M0yr“^ (Natta et al. 2006) to 3 x lO“®M0yr“^ 
(Greene & Lada 1996; Luhman & Rieke 1999). Transition 
discs like LkGalS and UX TauA, which also have low ac¬ 
cretion rates for their disc masses, have such large gaps 
(> 25 AU) that multiple giant planets, spread out over a 
range of radii, have been invoked to explain the large op¬ 
tically thin regions (Zhu et al. 2011). In contrast, T Tauri 
discs that have formed a single giant planet are expected to 
clear much smaller gaps, particularly if the planet is located 
at small disc radii, e.g., at a Jupiter-like distance of 5 AU. A 
high mass giant planet (~ 5Mj) located at such a distance is 
expected to greatly reduce the stellar accretion rate (Lubow 
et al. 1999). Given its large disc mass (~ O.IM0), it seems 
possible that DoAr25 may have formed such a giant planet. 
Interestingly, the accretion rate estimates in the literature 
for DoAr25 differ by more than a factor of 10. Extreme vari¬ 
ability (by such a large factor) in the stellar accretion rate 
is unusual (e.g., Nguyen et al. 2009) and may be consistent 
with the presence of a massive close companion (e.g., Basri 
et al. 1997). 

Other sources near Oph may share these properties. For 
example, the MIR colors of RXJ1633.9-2242 (source #32 
in Cieza et al. 2010) are consistent with the presence of a 
disk cavity ~ 8AU in radius (Orellana et al. 2012). The 
source has a fairly large disk mass (Ma — 0.01-0.02 M0) 
and a low estimated accretion rate based on its Ha width 
(M» = 1O“^°M0 yr“^). The low accretion rate could be con¬ 
firmed with a more accurate accretion rate diagnostic (e.g., 
IR hydrogen line luminosities or U-band excess). 

Cavities 5AU or 8AU in radius would have been diffi¬ 
cult to detect in the submillimeter studies of nearby young 
stars reported to date. The smallest dust cavity reported 
in the SMA study of Andrews et al. (2011) has a radius 
of 15 AU (WSB60). Higher angular resolution observations 
(e.g., with ALMA) of DoAr25 and the other low accretion 
rate systems would explore whether these systems are like 
the known Oph transition discs in having inner cavities. 

Thus we find that the combination of stellar accretion 
rates and disc masses can help to distinguish among pro¬ 
posed scenarios for the origin of transition discs. Toward this 
end, we have attempted to control for the known dependen¬ 
cies of M* and Ma on age (e.g., Hartmann et al. 1998), M* 
on stellar mass (e.g., Muzerolle et al. 2003), and Ma on the 
presence of stellar or sub-stellar companions (e.g., Artymow- 
icz & Lubow 1996; Lubow et al. 1999) hy (i) comparing the 
properties of transition discs and normal T Tauri stars in 
the same star-forming region, (ii) correcting for the known 
dependence of M* on M*, and (in) excluding known bina¬ 
ries, respectively. With this approach, we find that including 
disc masses in the demographic analysis helps to accentuate 
decrements in the stellar accretion rates of transition disc 
subpopulations relative to normal T Tauri stars. 

Our results complement previous studies that have 
taken different approaches in examining the accretion rates 


of transition discs in other star forming regions. Some of 
the earlier studies focused on sources older than those stud¬ 
ied here. For example, in the Tr 37 and NGG 7160 clusters 
in Cep OB2 (^ 4 Myr and ~ 12 Myr respectively), half 
the transition discs exhibit no measurable accretion, and 
the accretion rates of the other half are similar to those of 
classical T Tauri stars in the region (2 — 3 x 1O“®M0 yr“^; 
Sicilia-Aguilar et al. 2010). These results may differ from 
our result because the dominant physical mechanism that 
produces transition disc SEDs may differ between ages of 
1 Myr and 10 Myr; e.g., photoevaporation may dominate at 
late times. 

Several studies have compared the stellar accretion 
rates (but not disc masses) of transition discs and T Tauri 
stars in star forming regions similar in age to those studied 
here (< 3Myr). Similar in spirit to our results, some studies 
have reported evidence that transition discs have lower ac¬ 
cretion rates (Espaillat et al. 2012) or lower accretion lumi¬ 
nosities (Salyk et al. 2013) than normal T Tauri stars. Kim 
et al. (2013) found that transition discs in Orion A have 
accretion rates a factor of 10 lower, on average, than Tau¬ 
rus T Tauri stars of the same stellar mass. In contrast, Fang 
et al. (2013) found that accreting transition discs in L1641 
have a median accretion rate (from Ha) similar to that of 
the T Tauri stars in the region. However, their transition 
disc accretion rates determined from H/3 luminosity do have 
a lower median value (by a factor ~ 3) than T Tauri stars 
of the same stellar mass. H/3 luminosity may more reliably 
trace accretion rate than Ha, because Ha is more optically 
thick and it can be reduced in strength by wind absorption 
components (e.g., Muzerolle et al. 1998). 

The choice of comparison samples may also play a 
role in the the interpretation of accretion rates. Manara et 
al. (2014) compared the accretion rates of transition discs 
in a variety of star forming regions (Taurus, Ophiuchus, 
Ghamaeleon, Serpens) with the accretion rates of T Tauri 
stars in Lupus I and II. The transition discs appear biased 
toward lower M*, with a median M* ~ 3 times lower than 
the average rate for the T Tauri stars for the same M*. The 
use of Lupus T Tauri stars as the comparison sample may 
underestimate the M* decrement, because stellar accretion 
rate declines with age (Hartmann et al. 1998), and the Lu¬ 
pus I and II regions (~ 3 Myr; Alcala et al. 2014) may be 
older than most of the transition discs in the Manara et al. 
study (1-2 Myr for Tau, Oph, Gha; e.g., Furlan et al. 2009). 
The average M* for Lupus T Tauri stars (0.4-0.8 Mq) is 
1.3 X lO“®M0yr“^ (Alcala et al. 2014), whereas the mean 
M* for Taurus T Tauri stars is ~ 10“® Mq yr“^ (Hartmann 
et al. 1998). Thus, other studies in the literature appear to 
find results consistent with those reported here when we con¬ 
sider the accretion rate tracers and comparison populations 
used. 


5 SUMMARY AND CONCLUSIONS 

We find that the combination of stellar accretion rates and 
disc masses can aid in distinguishing among proposed ex¬ 
planations for the unusual SEDs of the population of young 
stellar objects known as transition discs. Determining how 
stellar accretion rate depends on a parameter such as SED 
shape can be a challenge because accretion rate is known to 



depend on other parameters such as age (e.g., Hartmann et 
al. 1998), stellar mass (e.g., Muzerolle et al. 2003), disc mass 
(Najita et al. 2007), and potentially the presence of stellar 
or sub-stellar companions (e.g., Artymowicz & Lubow 1996; 
Lubow et al. 1999). In our study, we have attempted to 
control for these known dependencies hy (i) comparing the 
properties of transition discs and normal T Tauri stars in 
the same star-forming region, (ii) correcting for the known 
dependence of M* on M*, (in) using available Md measure¬ 
ments, and (iv) excluding known binaries, respectively. 

Using this approach, we compared the distributions in 
the M*-Md plane of transition discs and single, normal T 
Tauri stars in Taurus and Ophiuchus. The transition discs 
we studied have prominent MIR dips in their Spitzer SEDs 
or inner cavities in submillimeter continuum images. The 
available data suggest some interesting trends: 

1. When combined with the Tau sample, the Oph sam¬ 
ple enhances the statistical signficance of the previously re¬ 
ported trend of M* vs. Md among single, normal T Tauri 
stars. 

2. The Tau and Oph transition discs we studied favor 
the lower right region of the M*-Md plane, i.e., they tend 
to have lower M* for a given Md and a higher Md (for a 
given M*) than normal T Tauri stars. This result is most 
consistent with the formation of objects massive enough to 
alter the accretion flow, i.e., single or multiple giant planets. 

3. Several T Tauri stars that are not known transi¬ 
tion disc systems also have low accretion rates for their 
disc masses. Such sources that also have high extinctions 
may have underestimated stellar accretion rates (LFAM3, 
WL14). The others may have formed (high mass) giant plan¬ 
ets that are massive enough to greatly reduce the stellar 
accretion rate, but the planets are located at orbital radii 
too small to have an observable impact on the submillime¬ 
ter continuum image or a distinctive impact on the SED 
(DoAr25, GY284). More detailed studies of the extinction 
toward these sources, their stellar accretion rates, and disc 
structure would help to distinguish between these possibili¬ 
ties. 
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Table 1. Ophiuchus Stellar Accretion Rates and Disk Masses 


Source 

(1) 

SpT 

(2) 

Mult 

(3) 

A* 

(4) 

logR* 

(5) 

logMf 

(6) 

logMf 

(7) 

logM^ 

(8) 

(9) 

Ref 

(10) 

Av 

(11) 

Type 

(12) 

Other Name 
(13) 

IRS2 

K3.5 

B 

1.66 

0.29 

0.16 

< -7.29 

< -3.58 



16.9 

N 

YLW19 

SR4 

K4.5 

S 

1.31 

0.28 

0.08 

-7.28 

-2.37 



2.2 

N 

AS206/IRS12 

GSS26 

K7 

S 

8.14 

0.76 

-0.09 

< -7.73 

-2.13 



21.9 

N 

- 

EL18 

Ml 

B 

2.14 

0.55 

-0.33 

< -8.01 

< -3.14 



9.9 

N 

GSS29 

DoAr24 

K5 

S 

1.23 

0.29 

0.02 

-8.88 

< -2.59 



3.5 

N 

GSS28/HBG638 

EL20 

MO 

S 

1.97 

0.50 

-0.27 

-6.82 

-2.15 

-1.54 

A09 

14.2 

N 

VSSGl 

GYll 

M6 

s 

0.002 

-0.78 

-1.13 

<-12.68 

< -2.48 



4.2 

N 

- 

DoAr25 

K5 

s 

1.52 

0.33 

0.04 

< -9.55 

-1.71 

-0.87 

A09 

3.4 

O 

WSB29/GY17 

LFAM3 

MO 

s 

0.31 

0.09 

-0.24 

-8.90 

-1.83 



16.1 

O 

GY21 

EL24 

K5 

s 

5.10 

0.60 

0.04 

-6.43 

-1.65 

-0.93 

A09 

8.7 

N 

WSB31 

VSSG27 

K6 

s 

1.83 

0.40 

-0.04 

-7.02 

< -2.82 



24.0 

N 

GY51 

WSB37 

M5 

B 

0.21 

0.19 

-0.69 

< -9.58 

< -2.42 



2.5 

N 

GY93 

IRS 17 

K3 

s 

1.14 

0.20 

0.17 

< -9.33 

< -2.58 



20.3 

N 

WL7/GY98/WLY2-17 

EL26 

MO 

s 

1.71 

0.47 

-0.28 

-8.04 

-2.93 



7.5 

N 

GSS37AB/GY110 

EL27 

MO 

s 

1.28 

0.40 

-0.28 

-7.10 

-1.57 

-0.84 

A09 

14.8 

N 

GSS39/GY116 

WSB38 

G3.5 

B 

5.46 

0.36 

0.07 

-8.64 

< -2.95 



7.0 

N 

VSS27AB/ROXsl6 

WL18 

K6.5 

s 

0.27 

0.00 

-0.10 

-7.84 

-2.45 



10.5 

N 

GY129 

WL14 

M4 

s 

0.14 

0.06 

-0.62 

<-10.20 

-2.28 



15.8 

O 

GY172 

WL16 

Al 

s 

250. 

0.79 

0.65 

-5.26 

< -3.93 



31.0 

N 

GY182/YLW5 

IRS29 

M4 

B 

0.63 

0.39 

-0.58 

< -9.90 

-3.18 



17.8 

N 

WL1AB/GY192 

WLIO 

K7 

s 

0.96 

0.29 

-0.12 

-8.63 

< -3.03 



11.5 

N 

GY211 

SR21 

F4 

s 

16.54 

0.49 

0.16 

< -8.51 

-2.40 

-2.22 

All 

6.9 

HT 

EL30/VSSG23 

WLll 

MO 

s 

0.10 

-0.15 

-0.27 

-10.33 

< -2.71 



11.5 

N 

GY229 

IRS34 

MO 

B 

1.07 

0.36 

-0.27 

< -8.41 

< -2.74 



24.7 

N 

GY239 

IRS39 

Ml.5 

s 

1.05 

0.41 

-0.38 

< -8.97 

-2.62 



18.4 

N 

WL4/GY247 

IRS42 

K7 

s 

7.29 

0.73 

-0.10 

< -7.57 

< -2.51 



25.0 

N 

GY252 

YLW16C 

Ml 

s 

1.39 

0.45 

-0.34 

-7.49 

-2.30 



22.5 

N 

GY262 

EL31 

M4 

B 

0.87 

0.46 

-0.58 

< -8.59 

< -3.02 



11.7 

N 

WL13/VSSG25/GY267 

EL32 

K6.5 

s 

1.75 

0.41 

-0.07 

< -8.65 

< -2.41 



21.1 

N 

IRS45/GY273 

EL33 

M3 

B 

4.88 

0.80 

-0.47 

< -7.70 

< -2.94 



24.1 

N 

IRS47/GY279 

GY284 

M3.2 

S 

0.12 

0.00 

-0.56 

<-10.31 

-1.62 



6.0 

O 

- 

ROXs25 

K7 

s 

3.07 

0.55 

-0.11 

-7.20 

-2.71 



11.6 

N 

GY292 

IRS49 

K5.5 

s 

2.16 

0.43 

0.22 

-8.30 

-2.70 



10.5 

N 

GY308 

WSB52 

K5 

s 

0.60 

0.27 

-0.34 

-8.26 

-2.15 

-2.15 

AlO 

5.0 

N 

GY314 

SR9 

K5 

B 

2.10 

0.40 

0.04 

-8.43 

< -3.12 



1.6 

N 

AS207/GY319/IRS52 

SRIO 

M2 

s 

0.33 

0.18 

-0.43 

-8.30 

< -2.49 



1.3 

N 

GY400/HBG265 

WSB60 

M4.5 

s 

0.20 

0.16 

-0.65 

-8.59 

-1.74 

-1.55 

All 

3.9 

H 

YLW58/B162816-243657 

SR20 

G7 

B 

8.71 

0.49 

0.07 

< -8.41 

< -3.30 



6.0 

N 

HBG643 

SR13 

M4 

B 

0.40 

0.29 

-0.59 

-8.46 

-2.09 

-1.92 

AlO 

0.0 

N 

HBG266 

AS205A 

K5 

S 

4.00 

0.54 

0.04 

-6.47 

-1.71 

-1.54 

AlO 

2.9 

N 

- 

AS209 

K5 

s 

1.50 

0.33 

0.04 

-7.29 

-1.76 

-1.60 

AlO 

0.9 

N 

- 

DoAr44 

K3 

s 

1.90 

0.31 

0.20 

-8.04 

-2.11 

-2.15 

All 

3.3 

HT 

ROXs44 

SR24S 

K2 

s 

4.40 

0.46 

0.30 

-7.15 

-1.88 

-1.35 

All 

7.0 

H 

- 

WaOph6 

K6 

s 

2.90 

0.50 

-0.04 

-6.64 

-2.02 

-1.11 

AlO 

3.6 

N 

- 


Column 2: Spectral types from McClure et al. (2010), Andrews Sz Williams (2007a, 2007b), Andrews et al. (2009, 2010). Column 3: 


Sources with a companion within l" are classified as binaries (Reipurth Sz Zinnecker 1993, Barsony et al. 2005; Ratzka et al. 2005). 
Columns 4 and 12: Stellar luminosity (in Lq) and Ay from McClure et al. (2010) and M. McClure (2014, private communication) 
except for WL16 (Ressler & Barsony 2003), WL14 (Luhman Sz Rieke 1999), GYll (Cushing et al. 2000), EL24, WSB52, SR13, AS205A, 
AS209, SR24S, WaOphO (Andrews et al. 2010). Columns 5 and 6: Derived stellar radius (in Rq) and stellar mass (in Mq). Column 
7: Stellar accretion rate derived from Natta et al. (2006) measurements using the stellar parameters from McClure et al. (2010) (see 
text for details). Column 8: Disk mass (in Mq) estimated from optically thin, isothermal approximation (see §2 for details). Column 9: 
Alternate disc mass estimate from a fit to the SED and resolved millimeter visibilities. Column 10: Reference for A09 = Andrews 

et al. (2009); AlO = Andrews et al. (2010); All = Andrews et al. (2011). Column 12: Source type, where N= normal T Tauri star; T= 
transition disc SED; H= submillimeter cavity; 0= outlier. 






Table 2. Taurus Stellar Accretion Rates and Disk Masses 


Source 

(1) 

SpT 

(2) 

logM* 

(3) 

logM* 

(4) 

log A4d 
(5) 

logMd.flt 

(6) 

Ref Md,fit 
(7) 

Type 

(8) 

AATau 

K7 

-0.10 

-8.48 

-2.04 

-1.52 

A07 

N 

BPTau 

K7 

-0.10 

-7.54 

-2.25 



N 

CITau 

K7 

-0.10 

-7.19 

-1.84 

-1.40 

A07 

N 

CWTau 

K3 

0.20 

-7.61 

-2.27 



N 

CXTau 

M2.5 

-0.44 

-8.97 

-2.77 



N 

CYTau 

M1.5 

-0.36 

-8.12 

-1.85 



N 

DETau 

Ml 

-0.41 

-7.59 

-2.39 



N 

DGTau 

K7 

-0.04 

-6.30 

-1.43 



N 

DHTauA 

Ml 

-0.32 

-8.30 

-2.54 

-2.52 

A07 

N 

DKTau 

MO 

-0.15 

-7.42 

-2.69 



N 

DLTau 

K7 

-0.09 

-6.79 

-1.61 

-1.00 

A07 

N 

DMTau 

Ml 

-0.33 

-7.95 

-1.72 

-1.40 

All 

T 

DNTau 

MO 

-0.22 

-8.46 

-1.93 

-1.22 

A07 

N 

DOTau 

MO 

-0.33 

-6.85 

-1.88 



N 

DRTau 

K7 

0.08 

-6.50 

-1.89 

-2.00 

A07 

N 

DSTau 

K5 

0.02 

-7.89 

-2.62 



N 

FMTau 

MO 

-0.22 

-8.45 

-2.63 



N 

FYTau 

K7 

0.07 

-7.41 

-2.80 



N 

FZTau 

MO 

-0.35 

-7.32 

-2.77 



N 

GITau 

K7 

-0.10 

-7.69 

-2.81 



N 

GKTau 

K7 

-0.11 

-8.19 

-3.44 



N 

GMAur 

K3 

0.13 

-8.02 

-1.66 

-1.15 

All 

T 

GOTau 

MO 

-0.20 

-7.93 

-1.98 

-0.74 

A07 

N 

HKTau 

M0.5 

-0.27 

-7.27 

-2.27 



N 

HNTau 

K5 

-0.04 

-8.41 

-2.64 



N 

HOTau 

M0.5 

-0.25 

-8.86 

-2.34 



N 

IPTau 

MO 

-0.22 

-9.10 

-2.83 



N 

IQTau 

M0.5 

-0.27 

-7.55 

-2.06 



N 

LkCal5 

K5 

0.02 

-8.87 

-1.74 

-1.26 

All 

T 

RWAurA 

K3 

0.18 

-7.12 

-2.52 



N 

RYTau 

K1 

0.44 

-7.11 

-1.92 

-1.70 

A07 

N 

UXTauA 

K2 

0.20 

-9.00 

-2.22 

-2.15 

All 

T 

V836Tau 

K7 

-0.08 

-8.98 

-2.33 



N 

Column 3: 

Stellar 

mass (in 

Mq). 

Column 4: 

Stellar accretion rate (in 

Mq yr' 


Column 5, 6: Disk mass (in Mq). Column 7; A07 = Andrews Williams (2007b); All 
= Andrews et al. (2011). Column 8: Source type, where N= normal T Tauri star; T= 
transition disc. 



